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摘 要 : 为 了 实现 基于 单独 光学 遥感 数据 对 科尔沁 沙丘 草 甸 相 间 地 区 土地 利用 与 覆 被 (LULC) 类 型 的 识别 ,选用 


2018 年 64 景 Sentinel-2 影 像 , 结 合影 像 分 割 技术 ,利用 植被 物候 信息 和 生境 特征 ,建立 了 基于 群落 水 平 的 LULC 决 


策 树 识 别 规则 ,总 体 分 类 精度 为 0.91,Kappa 系 数 为 0.89。 分 类 结果 显示 :研究 区 旱地 分 布 面积 最 大 , 占 33.79% HE 


木 群落 次 之 , 占 25.03%, 高 多 样 性 半 灌 木 群 落 和 乔木 林 相近 ,分 别 为 14.54% 和 10% , 低 多 样 性 半 灌 木 群 落 . 草 多 地 和 


流动 沙 地 分 别 占 5% 左 右 , 剩 余 类 型 的 总 占 比 小 于 5%。 该 方法 不 仅 可 以 准确 反映 研究 区 覆 被 类 型 的 空间 分 布 情 
况 , 还 能 给 出 不 同 覆 被 类 型 的 生长 发 育 状况 ,可 为 该 区 域 物质 循环 研究 提供 基础 数据 ,同时 为 该 区 域 历史 LULC 识 


别提 供 阔 值 参考 。 
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科尔沁 沙 地 是 中 国 四 大 沙 地 中 最 大 的 沙 地 ,其 
生态 系统 脆弱 ,是 华北 地 区 沙尘暴 的 主要 来 源 "?”。 
土地 利用 与 履 被 (land use/land cover, LULC ) 监 测 是 
该 地 区 风沙 化 趋势 预测 的 基础 。Duan 等 中 利用 
Landsat 系列 数据 研究 了 当地 土地 风沙 化 的 时 空 动 
态 , 指 出 2010 年 的 沙化 土地 占 监测 面积 (1.25x105 
km ) 的 24%, 而 71% 为 轻 度 和 中 度 沙 化 土地 ,植被 覆 
盖 度 达到 30% 以 上 ,主要 分 布 在 科尔沁 沙 地 的 中 部 
和 东部 。 位 于 科尔沁 沙 地 东南 边缘 , 科 左 后 旗 的 沙 
化 土地 占 监测 面积 (11632 km2) 的 50% ,其 中 78% 为 
轻 度 和 中 度 沙 化 土地 ,可 见 当 地 的 生态 修复 措施 起 
到 了 显著 成 效 “。 男 外 ,有 研究 表明 沙 地 的 LULC 和 
水 量 平衡 具有 复杂 关系 ,例如 ,在 科尔沁 沙 地 南 
部 的 张 古 台 , 由 于 阔 叶 林 及 农田 种 植 面积 的 增加 ， 
导致 当地 的 地 下 水 位 下 降 ,从 而 使 原始 引进 的 樟 子 
松 人 工 林 出 现 枯 奢 ,出 现 大 面积 死亡 的 现象 。 
此 ,识别 科尔沁 沙 地 的 LULC 对 当地 的 土地 管理 调 
配 及 可 持续 性 发 展 研究 具有 重要 意义 。 

随 着 遥感 数据 源 的 不 断 丰 富 以 及 遥感 应 用 和 
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计算 机 技术 的 快速 发 展 ,利用 遥感 数据 获取 土地 利 
用 与 覆 被 信息 已 成 为 识别 LULC 的 主要 手段 5 。 前 
期 关于 科尔沁 沙 地 LULC 的 研究 多 数 是 基于 植被 覆 
善 度 进 行 的 分 类 。 随 着 当地 沙化 土地 的 不 断 修 
复 ,我 们 研究 的 重心 应 该 转移 到 植被 覆盖 类 型 的 精 
细 划 分 以 及 植被 多 样 性 监测 上 。He 等 "指出 监 
测 外 来 物种 对 生境 结构 和 条 件 改造 以 及 关键 物种 
发 生 和 传播 变化 的 预警 信号 对 于 生态 系统 的 管理 
至 关 重 要 。 男 外 ,Rocchini 等 ”通过 对 比 IKONOS 
和 Landsat 数 据 反 演 植 被 种 类 及 多 样 性 时 的 精度 , 指 
出 对 于 植被 多 样 性 的 反 演 ,光谱 分 状 率 比 空间 分 辨 
率 更 重要 。Conese FI Maselli "指出 在 复杂 地 貌 的 
植被 分 类 过 程 中 可 以 通过 多 时 相 数 据 ,借助 植被 物 
候 信息 从 而 提高 反 演 分 类 精度 。Sentinel-2 卫星 是 
欧 空 局 哥 白 尼 计 划 中 发 射 的 第 2 组 卫星 , 它 不 仅 比 
Landsat 卫星 产品 的 空间 分 辩 率 高 ,而且 具 有 包括 可 
见 光 、 近 红外 (NIR) 植被 红 边 和 短波 红外 (SWIR ) 
区 的 13 个 光谱 波段 。 男 外 ,还 具有 重 访 周期 短 达 五 
天 的 高 时 间 分 辨 率 ,因此 ,使 用 多 时 相 Sentinel- 
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2 卫星 数据 有 望 实现 群落 水 平 上 的 LULC 精准 识别 。 

本 文 以 科尔沁 东南 边缘 的 沙丘 草 负 相间 地 区 
为 例 , 基 于 2018 年 植被 生长 季 64 景 Sentinel-2 遥感 
数据 和 地 面 调查 数据 ,利用 植被 物候 信息 和 生境 特 
征 ,结合 基于 面向 对 象 影像 分 割 技术 的 多 分 辨 率 分 
制 (Multi-Resolution Segmenation, MRS ) 算 法 ,建立 基 
于 群落 水 平 的 LULC 决 策 树 识 别 规则 ,以 求 给 出 能 
反应 下 垫 面 植被 生长 发 育 信息 的 研究 区 LULC 空间 
分 布 格局 ,为 后 期 该 地 区 的 生态 系统 健康 评价 提供 
参考 数据 。 


1 研究 区 概况 


人 研究 区 位 于 科尔沁 沙 地 东南 缘 , 是 西 辽 河流 域 
中 的 一 个 财 流 区 域 。 隶 属于 内 蒙古 自治 区 通辽 市 
科尔沁 左翼 后 旗 , 为 典型 的 荒漠 化 的 沙丘 草 旬 相间 
地 区 , HOF AA HE 122°10’~123°10'E, 42°50'~43°46'N, 
海拔 高 度 133~267 m ,总 面积 约 $787.32 km’, , 属 半 干 
旱 - 半 湿 润 大 陆 性 季风 气候 。 人 研究 区 冬季 主要 受到 
蒙古 地 区 冷 高 压 的 控制 ,多 刮 偏 北 风 或 偏 西风 ,在 
夏季 主要 受 大 陆 低 气压 以 及 副热带 高 压 所 控制 ,多 
以 偏 南 和 西南 风 为 主 。 地 表 水 和 地 下 水 资源 较 丰 
富 ,湖泊 沼泽 ,水 泡 密 布 ( 图 1) ,整个 地 表 由 多 个 相 
似 小 单元 随机 组 合 形 成 微波 状 起 伏 平原 ,每 个 小 单 
元 是 由 包围 着 湖泊 的 草 匈 斑 块 和 其 周围 的 沙丘 斑 
块 组 成 的 中 。 本 文选 择 研 究 区 中 的 一 个 典型 小 单 
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元 作为 LULC 识 别 决策 树 规则 建立 的 试验 区 。 
2 数据 来 源 


2.1 Sentinel-2 遥 感 影像 

本 研究 使 用 Sentinel-2 的 level 1C 级 影像 产品 
(经 过 正 射 校正 和 亚 像 元 级 几何 精 校 正 的 大 气 表 观 
反射 率 ) 作 为 原始 分 析 数 据 ,影像 通过 欧洲 空间 局 
Px} 44 Chttps://scihub.copernicus.eu ) #7 38 3k A, FETT 
括 2018 年 植被 生长 季 16 次 过 境 的 64 景 产品 。 图 2 
列 出 了 过 境 的 具体 时 间 和 试验 区 的 气候 状况 。 使 
用 Sentinel Application Platform (SNAP) 平 台 和 
Sen2cor2.4 插 件 " 对 level 1C 级 产品 进行 辐射 定 标 
和 大 气 校正 生成 L2A 级 数据 。 采 用 最 邻近 内 插 法 
将 L2A 级 数据 的 10 个 波段 重 采样 生成 空间 分 辨 率 
为 10 m 的 产品 ,10 个 波段 分 别 是 2~4,8(10 m) 、5~ 
7,8A,11 和 12(20 m)。 然 后 利用 SNAP 软 件 对 各 个 
时 间 的 每 个 波段 的 4 景 影像 进行 馈 般 ,最 后 利用 
ENVI 软 件 将 产品 裁剪 到 研究 区 。 
2.2 地 面 调 查 数 据 

2018 年 7 一 9 月 ,根据 沙丘 草 僻 相间 分 布 的 地 貌 
类 型 ,在 试验 区 布设 了 54 个 调查 点 。 在 调查 时 ,对 
于 生长 分 布 整齐 的 农田 和 乔木 林 只 做 踩点 调查 ,对 
于 生长 分 布 不 规则 的 其 他 植被 做 精细 的 调查 。 在 
每 个 精细 调查 点 先 布设 一 个 10 mx10 m 的 灌木 . 半 
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图 1 研究 区 及 试验 区 位 置 示意 图 


Fig.1 The geographical location and the relief map of the study area and test area 
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图 2 2018 年 试验 区 气象 条 件 及 遥感 卫 片 过 境 时 间 


Fig. 2 Meteorological conditions and transit date of satellite in the test area of 2018 


灌木 调查 大 样 方 ,然后 在 大 样 方 的 4 角 和 中 心 共计 
布设 5 个 1 mxl m 的 草本 调查 小 样 方 。 分 别 调查 样 
方 内 植被 的 种 类 高度、 密度, 对 于 无 性 繁殖 物种 ， 
其 分 株 被 视 为 物种 密度 的 基本 单位 ”。 最 后 用 收 
割 法 分 物种 获取 样 方 内 每 一 种 植物 的 3 个 标准 株 ， 
装 和 人 密封 的 塑料 袋 带 回 实验 室 ,用 烘箱 60 CTH 
48h ,测定 其 干 质量 , 乘 以 密度 得 到 该 物种 的 地 上 生 
物 量 。 利 用 获得 的 生物 量 .密度 ,高度 计算 每 个 调 
查 点 的 Shannon 多样 性 指数 '”。 


3 研究 方法 


3.1 研究 区 LULC 分 类 

参照 前 期 对 试验 区 进行 的 植被 分 类 研究 ”2 ， 
结合 2018 年 地 面 调查 数据 ,研究 区 LULC 被 划分 为 
10 类 ,分 别 是 :(1) 居民 用 地 ,表示 房屋 分 布 聚集 的 
区 域 ;(2) 水 域 ;(3) 旱地 ,主要 作物 是 玉米 ;(4) 水 
田 ,主要 作物 是 水 稳 ; (5) 乔木 林 , 代表 树 种 是 杨 树 ; 
(6) 灌木 群落 ,优势 物 种 是 小 叶 锦 鸡 儿 (Caragana 


microphylla) ;(7) 低 多 样 性 半 灌 木 群落 ,优势 物种 是 
E BEIR S (Artemisia halodendron) ; (8) 高 多 样 性 半 
灌木 群落 ,优势 物种 有 冷 葛 (4rtemisia frigida) ER 
W (Ephedra sinica) 和 冰 草 (Agropyron cristatum ) ; (9) 
草 馈 地 ,优势 物种 有 蔡 草 (Carex duriuscula) 和 芦苇 
(Phragmites communis ); (10) 流动 沙 地 ,包括 无 植被 
图 盖 和 有 覆盖 少量 植被 两 种 情况 , 当 有 植被 覆盖 时 ， 
优势 物种 是 沙 米 (hgriophyllum squarrosum ) 和 差 巴 
叮 草 ,群落 中 植被 种 类 少 ,生物 量 和 多 样 性 低 。 不 
同 植被 群落 类 型 的 组 成 及 结构 特征 信息 见 表 1。 
3.2 研究 方法 概述 

基于 群落 水 平 研究 区 LULC 识别 步骤 如 下 :第 
一 步 基 于 试验 区 实测 数据 构建 LULC 决 策 树 识别 规 
则 。 首 先 探寻 影像 最 优 分 割 方式 ,然后 导入 训练 样 
本 生成 对 象 样本 ,分 析 不 同 LULC 类 型 的 归 一 化 植 
被 指数 (Normalized Difference Vegetation Index, ND- 
V1) 和 各 个 波段 反射 率 的 时 序 变 化 规律 ,建立 识别 
规则 ,并 根据 验证 样本 调试 识别 规则 使 解 译 结果 最 
接近 真实 情况 。 其 中 训练 样本 为 2018 年 的 地 面 实 


表 1 土地 利用 / 覆 被 类 型 的 实测 数据 统计 结果 


Tab.1 Statistical results of field data on land use/cover types 


土地 利用 / 歼 被 类 型 优势 物种 调查 点 /个 ”植被 种 类 (类 ) 地 上 单位 面积 生物 量 /(g'm”) ”Shannon 多样 性 指数 (H) 验证 样本 
草包 地 FE aE 6 9.17+4.95 438.65+340.23 1.51+0.68 125 
流动 沙 地 Wok FE 5 6.33+1.6 107.35+58.08 1.09+0.31 55 
REREN HERS 4 9.25+1.09 122.55+77.17 1.66+0.10 75 
ABEREEK YES JRE RE 7 13.67+2.92 229.32+155.94 1.99+0.13 110 
灌木 群落 小 叶 锦 鸡 儿 16 13.88+4.57 1292.5+1367.16 1.81+0.25 130 
乔木 林 杨 树 6 - - 70 
旱地 玉米 7 一 一 一 185 
水 | 水 稻 3 一 一 一 25 


注 :- 表 示 无 测量 数据 。 
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测 数据 ,验证 样本 是 2015 年 采集 的 31 个 50 mx50 m 
的 调查 点 ,具体 样本 数量 分 配 见 表 1。 对 于 水 域 
和 有 居民 用 地 特征 明显 的 地 物 ,采用 人 工 勾 绘 的 方法 
进行 验证 样本 添加 。 第 二 步 应 用 建立 的 识别 规则 
实现 研究 区 LULC 解 译 ,并 借助 Google earth 对 结果 
进行 粗 评 。 
3.2.1 多 尺度 影像 分 割 ” 解 译 过程 中 的 影像 分 割 通 
过 多 分 状 率 分 割 (MRS) 算 法 实现 ”。 该 算法 具有 3 
个 用 户 定义 参数 :尺度 、 形 状 和 紧凑 度 , 尺 度 参 数 是 
判断 所 考察 像素 是 否 合 并 到 临近 影像 对 象 的 标准 ， 
形状 参数 的 值 等 于 1 减 去 颜色 参数 , 紧凑 度 参 数 的 
值 等 于 1 减 去 光滑 度 参 数 ”。 通 过 不 同 分 割 参数 产 
生 的 不 同 对 象 层 之 间 存 在 属性 继承 关系 ,确定 好 分 
割 参数 对 于 LULC 的 正确 识别 很 重要 “””。 首 先 通 
过 目 视 法 ,选取 居民 用 地 、 水 域 和 植被 3 种 差异 显著 
的 地 物 挑选 用 于 影像 分 割 的 数据 源 ,并 采用 ESP2 
(Estimating the Scale Parameter 2) fx (ii K E TA” 
结合 人 工 调试 参数 对 比 的 方法 获得 最 优 分 割 参数 。 
3.2.2 生长 季 时 序 NDVI 数 据 处 理 NDVI 指数 是 目 
前 应 用 最 广泛 的 植被 指数 ,计算 公式 如 下 : 

NDVI = nr rep (1) 

R 


NIR + Reep 
式 中 :Rim 为 红 光 波段 的 反射 率 ; Rn 为 近 红 外 波段 
的 反射 率 ,分 别 为 Sentinel-2 数 据 的 波段 4 和 8。 
首先 通过 对 计算 生成 的 初始 NDVI 时 序数 据 进 
行 线性 插值 ,获得 4 月 4 日 至 10 月 31 日 每 间隔 15 d 
的 NDVI 时 序数 据 , 例 如 ,NDVLnmisa=NDVInsa+(ND- 
VLnan-NDVILnsn'(1O015)。 为 进一步 去 除 NDVI 时 
序数 据 中 的 噪声 ,使 其 更 为 符合 植被 生长 发 育 
过 程 ,采用 双重 逻辑 困 数 滤波 (Double Logistic , 
DD) 法 天 对 NDVI 时 序数 据 进行 平滑 处 理 ( 图 3), 采 
用 seasonal amplitude 方 法 模拟 植被 生长 过 程 ,将 
Start time 和 End time 的 参数 均 设 置 为 0.5 ,获得 植被 
NDVI 季节 最 大 值 (NDVI) 和 生长 季 开 始 时 间 
(Time.) ,以 上 操作 过 程 均 在 Timesat 3.1 软件 中 
完成 。 
3.2.3 非 植被 识别 “对 水 域 的 识别 使 用 归 一 化 水 体 
6 2#U (Normalized Difference Water Index,NDWT) , 计 


NDWI= Roacen 7 Rue (2) 
GREEN T RNIR 
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0.8 


0.6 


一 上 一 原始 数据 
去 品 数 据 
4  TIMEs 

= NDVImax 


0.2 


0 1 
04-04 05-04 06-03 07-03 08-02 09-01 10-01 10-31 
日 期 /月 -日 
图 3 2018 年 单个 像 元 乔木 林 时 序 NDVI 双 逻辑 
斯 带 滤波 处 理 过 程 
Fig. 3 The processed NDVI curves of single 


pixel about forest by Double Logistic filter 


TAP : Rore MARIE BE HY NZE 5 Run AUT ZT INI BS 
的 反射 率 , 分 别 为 Sentinel-2 数 据 的 波段 3 和 8。 

流动 沙 地 的 划分 使 用 根据 像 元 二 分 模型 推 
导出 的 植被 覆盖 指数 (Vegetation Coverage Index， 
VCI) ,计算 公式 如 下 : 


-NDVI -NDVI 
ar NDVI, a -NDVI oi (3) 


式 中 :NDVI。 为 裸 地 对 应 的 像 元 数值 ;NDVI.. 为 完 
全 被 植被 覆盖 时 的 像 元 数值 ,在 本 研究 中 ,NDVI 设 
置 为 NDVIL ,根据 地 面 调查 数据 ,将 NDVIL 设 置 为 
0.25 ,NDVI 设置 为 0.8。 

由 于 居民 用 地 中 包含 NDVI 较 低 的 裸 地 和 房 
屋 , 以 及 NDVI 较 高 的 绿化 和 院内 种 植 土地 ,所 以 通 
过 NDVI 值 及 波段 反射 率 的 时 序 变 化 很 难 划分 。 首 
先 基 于 8 月 10 日 的 真 彩 图 ,通过 房屋 的 颜色 选择 样 
本 , 即 选择 颜色 条 件 满足 Rz200,G<150,B<150 的 对 
象 为 样本 ,然后 在 影像 分 割 的 基础 上 ,通过 基于 样 
本 分 类 功能 划分 出 居民 用 地 ,划分 流程 如 图 4 所 示 。 


4 结果 与 分 析 


4.1 LULC 决 策 树 识别 规则 

图 Sa 为 生长 季 7 种 覆 被 类 型 的 NDVI 变 化 曲 
线 ,与 植被 的 物候 信息 相关 。 可 以 看 出 ,由 于 乔木 
林 和 草 多 地 的 返青 时 间 比 其 他 植被 早 , 从 5 月 初 到 6 
月 初 ,他 们 的 NDVI 均 值 高 于 其 他 植被 。 由 于 水 田 
和 旱地 的 返青 时 间 和 快速 生长 期 比 其 他 植被 晚 ,从 
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影像 分 割 


图 4 


基于 颜色 筛选 样本 


基于 样本 分 类 


导 民 用 地 遥感 解 译 步骤 示意 网 


Fig.4 Schematic diagram of village remote sensing interpretation steps 


1.0) (a) 7 种 覆 被 类 型 的 NDVI 平 均值 及 
最 大 值 最 小 值 时 序 变化 


NDVI 


04-04 05-04 06-03 07-03 08-02 09-01 10-01 10-31 
日 期 /月 -日 


一 之 一 水 田 


一 一 乔木 林 ”一 一 旱地 


0.6 [ (b) 短波 红外 波段 平均 反射 率 的 时 序 变化 


04-04 05-04 06-03 07-03 08-02 09-01 10-01 10-31 
日 期 /月 -日 


一 一 高 多 样 性 半 灌 木 群落 


一 一 低 多 样 性 半 灌 木 群 落 。 一 一 灌木 群落 ”一 一 草 甸 地 


图 5 各 和 覆 被 类 型 的 遥感 反 演 信息 特征 


Fig.5 Features of remote sensing inversion information of various types of cover 


~ 


月 下 旬 到 10 月 初 , 他 们 的 NDVI 均 值 高 于 其 他 植 
被 。 由 于 灌木 的 返青 时 间 比 半 灌 木 早 且 长 势 更 成 
成 ,在 研究 期 间 ,灌木 群落 的 NDVI 均 值 大 于 半 灌 木 
群落 。 图 5b 为 不 同 覆 被 类 型 对 短波 红外 波段 (B11) 
反射 率 的 时 间 变 化 序列 ,与 植被 下 垫 面 的 水 分 条 件 
有 关 。 在 研究 期 间 ,水 田 的 平均 反射 率 值 最 低 , 草 
旬 地 次 之 , 低 多 样 性 半 治 木 群 落 最 高 ,其 他 履 被 类 
型 的 反射 率 比 较 接近 。 虽 然 各 种 覆 被 类 型 的 遥感 
特征 在 平均 水 平 上 差异 较为 显著 ,但 是 由 于 相同 覆 
被 类 型 之 间 存 在 一 定 的 发 育 差异 ,因此 ,需要 寻找 
合适 的 划分 规则 和 阅 值 使 解 译 结 果 最 接近 真实 
情况 。 


综 上 所 述 ,经 过 反复 调试 对 比 ,建立 的 研究 区 
LULC 决 策 树 识别 规则 如 图 6 所 示 。 首 先 划 分 居民 
用 地 ,基于 8 月 10 日 影像 的 10 个 波段 进行 第 一 次 分 
割 , 分 制 尺 度 、 形 状 参数 、 光 滑 度 分 别 设置 为 108、0.9 
和 0.1, 基 于 导入 的 房屋 样本 进行 分 类 。 然 后 划分 水 
域 ,基于 5 月 4 日 影像 的 10 个 波段 进行 第 二 次 分 割 |， 
AY HIRE .形状 参数 .光滑 度 分 别 设 置 为 83.0.4 和 
0.5 ,为 了 避免 水 中 植被 对 NDWI 指 数 的 影响 ,以 及 
考虑 到 研究 区 4 月 19 日 影像 部 分 有 云 覆 盖 的 情况 ， 
将 4 月 19 日 和 5 月 4 日 两 期 产品 的 NDWI 指 数 的 最 
大 值 大 于 0 的 对 象 定义 为 水 域 。 接 着 对 研究 区 的 覆 
被 类 型 进行 详细 识别 ,基于 8 月 10 日 影像 的 10 个 波 
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第 一 次 数据 源 : 8 月 10 日 的 10 个 波段 
分 割 分 割 尺度 : 108 
形状 参数 : 0.9 光滑 度 : 0.1 


第 二 次 数据 源 : 5 月 4 日 的 10 个 波段 
Ker 分 割 尺 度 : 83 
形状 参数 : 0.4 光滑 度 : 0.5 


第 三 次 数据 源 : 8 月 10 日 的 10 个 波段 
Aya 分 割 尺度 : 52 
形状 参数 : 0.4 光滑 度 : 0.5 


WV oe 
未 分 类 是 否 
,i 
未 分 类 
as 一 
未 分 类 

et T 是 


图 6 研究 区 土地 利 月 


低 多 样 性 半 高 多 样 性 半 
灌木 群落 灌木 群落 
且 履 被 类 型 决策 树 分 类 


Fig. 6 Classification decision tree of land use/cover type in the study area 


段 进行 第 三 次 分 割 ,分 割 尺度 .形状 参数 .紧凑 度 分 
别 设置 为 52、0.4 和 0.5。 首 先 提 取 流 动 沙 地 ,将 VCI 
小 于 10% 的 对 象 定义 为 流动 沙 地 。 其 次 提取 农田 ， 
将 返青 时 间 大 于 国 值 6(4 月 4 日 为 1, 每 隔 15 d 数 值 
加 1), 上 且 7 月 上 名 的 NDVI 增 长 率 大 于 阔 值 0.1 的 对 
eA AL NDV Inn ACF REL 0.5 的 对 象 定义 为 农田 。 
接着 提取 草 多 地 ,因为 研究 区 草 向 地 表层 含水 量 比 
固定 、 半 固定 沙丘 \ 流 动 沙 丘 普遍 偏 大 的 特点 "3, 将 
对 植被 和 土壤 水 分 敏感 的 短波 红外 波段 的 反射 率 
小 于 某 个 阔 值 的 对 象 定义 草 甸 地 ,采用 5 月 24 日 影 
像 的 B11 数 据 , 赣 值 为 0.3。 接 着 提取 乔木 林 ,将 返 
青 时 间 小 于 特定 阔 值 5 的 对 象 定义 为 乔木 林 。 最 后 
JEP SE JE BE pe FBC RS) FE OH LU R EE ACRES 
和 灌木 群落 ,定义 在 6 月 3 日 NDVI 的 模糊 区 间 从 小 
增长 到 大 是 半 灌 木 群 落 到 灌木 群落 ,模糊 区 间 设 定 


为 0.25~0.3。 

使 用 对 水 汽 变化 比较 敏感 的 Bl11 将 农田 划分 为 
旱地 和 水 田 ,与 草 多 地 识别 使 用 相同 的 数据 , 阔 值 
为 0.2。 因 为 分 布 高 多 样 性 半 灌 木 群 落 的 沙丘 恢复 
状态 比分 布 低 多 样 性 半 灌 木 群 落 的 好 ,土壤 的 持 水 
PERE IEC ,将 对 B11 反 射 率 大 于 某 个 国 值 的 对 
象 定 义 为 低 多 样 性 羊 灌木 群落 ,否则 是 高 多 样 性 
半 灌 木 群落 ,采用 的 是 8 月 2 日 的 数据 , 闪 值 设 定 
为 0.4。 

4.2 LULC 识别 结果 

研究 区 的 LULC 空间 分 布 如 图 7 所 示 。 本 人 研究 
的 总 解 译 面 积 为 5558.91 km , 其 中 旱地 的 分 布 面 
积 最 大 , 占 总 面积 的 33.79%。 灌 木 群 落 次 之 , 占 
25.03% , 主要 分 布 在 研究 区 的 西北 部 分 。 高 多 样 
性 半 灌 木 群落 和 乔木 林 的 分 布 面积 相近 ,分 别 占 总 
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Fig. 7 Interpretation results of land use/cover 


types in the study area 


面积 的 14.54% 和 10% 。 剩 余 LULC 的 分 布 面积 相 
对 较 小 , 占 总 面积 的 百分比 由 小 到 大 为 :水 田 
0.51% ,水 域 0.94% , 居民 用 地 2.91% ,流动 沙 地 
4.79% , WK & PE PE AE HE OAR BE YR 5.66% , 草 匈 地 
5.67% ,其 中 高 多 样 性 半 灌 木 群落 也 主要 分 布 在 研 
究 区 的 西北 部 分 , 低 多 样 性 半 灌 木 群 落 主 要 围绕 
流动 沙 地 分 布 。 


4.3 分 类 结果 的 评价 

为 了 验证 构建 分 类 规则 的 精度 和 应 用 潜力 , 首 
先 根据 验证 样本 对 试验 区 的 分 类 结果 进行 精度 评 
价 , 选 用 的 评价 指标 有 Kappa 系数、 总 体 分 类 精度 、 
用 户 精度 和 制图 精度 ,然后 结合 Coogle earth 数据 对 
研究 区 的 分 类 结果 进行 粗 评价 。 其 中 用 户 精 度 表 
示 在 被 分 为 某 类 的 所 有 样本 中 ,其 真实 样本 所 占 的 
比例 。 制 图 精度 表示 某 类 的 所 有 真实 样本 被 正确 
划分 的 比例 。 总 体 分 类 精度 表示 所 有 样本 中 被 正 
确 分 类 的 像 元 总 数 除 以 总 像 元 数 。Kappa 系 数 能 够 
结合 所 有 因素 综合 评价 分 类 的 结果 ,更 加 准确 的 反 
映 分 类 的 整体 精度 。 

从 表 2 中 可 见 , 训 练 样本 的 总 体 分 类 精度 为 
0.92 ,用 户 精 度 除 乔木 林 和 灌木 群落 以 外 均 在 0.90 
以 上 ,制图 精度 除 高 多 样 性 半 灌 木 群 落 以 外 都 在 
0.80 以 上 ,可 能 是 由 于 部 分 乔木 林 的 生长 状态 不 符 
合 正 常 物候, 例如 部 分 人 工种 植 的 杨 树 林 在 缺 水 状 
态 下 发 育 不 良 ,因此 可 能 被 解 译 为 灌木 群落 。 其 中 
水 域 . 居 民用 地 和 水 田 的 用 户 精 度 和 制图 精度 均 为 
1, 经 计算 Kappa 系数 为 0.89, 说 明 对 试验 区 的 分 类 
结果 是 比较 接近 真实 状况 的 。 另 外 ,对 于 人 研究 区 分 
类 结果 ,借助 Google earth 工具 ,随机 选择 研究 区 内 
的 3 个 位 置 ,对 范围 内 能 够 识别 的 水 域 流动 沙 地 、 
居民 用 地 .乔木 林 和 灌木 群落 进行 模糊 验证 ,同样 
达到 了 几乎 完全 一 致 的 评价 水 平 。 


表 2 试验 区 分 类 精度 


Tab. 2 Classification accuracy of study area 


WE FHE 高 多 样 性 济 云 

FAN pa yanma MAng TME DE AN pa 水城。 用 地 Ma 
乔木 林 60 0 0 0 25 0 0 0 0 0 0.71 
灌木 群落 10 130 0 28 0 10 0 0 0 0 0.73 
低 多 样 性 半 灌 木 群 落 0 0 65 0 0 0 0 0 0 0 1 
高 多 样 性 半 灌 木 群 落 0 0 5 82 0 0 0 0 0 0 0.94 
草包 地 0 0 0 0 100 0 0 0 0 0 1 
旱地 0 0 0 0 0 175 0 0 0 0 1 
水 田 0 0 0 0 0 0 25 0 0 0 1 
流动 沙 地 0 0 5 0 0 0 0 55 0 0 0.92 
水 域 0 0 0 0 0 0 0 0 50 0 1 
居民 用 地 0 0 0 0 0 0 0 0 0 75 1 
制图 精度 0.86 1 0.87 0.75 0.8 0.95 1 1 1 1 
总 体 分 类 精度 0.91 
Kappa 系 数 0.89 
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5 结论 


采用 时 间 分 辩 率 .空间 分 辨 率 和 光谱 分 辨 率 均 
较 高 的 Sentinel-2 数 据 , 根 据 植被 物候 信息 以 及 生 
境 条 件 对 研究 区 植被 进行 识别 。 依 据 返 青 时 间 , 覆 
被 类 型 可 划分 为 早 .中 、 晚 3 大 类 ,分 别 是 返青 最 早 
的 乔木 林 和 草 多 地 ,返青 最 晚 的 农田 以 及 中 间 返 青 
的 其 他 植被 群落 。 另 外 ,根据 植被 生长 环境 的 水 分 
ARE ,可 以 将 履 被 类 型 划分 为 种 植 在 水 中 的 水 稻 ， 
分 布 在 水 域 周 围 水 分 条 件 相 对 较 好 的 植被 ,有 玉米 
地 与 草 旬 群落 ,以 及 分 布 在 沙 地 水 分 条 件 相对 差 的 
植被 。 此 外 ,还 根据 沙 地 植被 多 样 性 高 低 与 土壤 表 
层 水 分 条 件 的 关系 ,对 半 灌 木 群 落 进 行 了 划分 。 根 
据 上 述 植被 特征 并 借助 多 分 辨 率 分 割 算法 建立 决 
策 树 识别 规则 ,总 体 分 类 精度 为 0.91,Kappa 系数 为 
0.89, 不 仅 达 到 了 人 研究 区 土地 利用 与 覆盖 类 型 精细 
划分 的 目的 ,而 且 能 够 反映 不 同 下 垫 面 的 物候 信息 
和 生境 条 件 ,为 以 后 研究 区 的 生态 系统 健康 评价 提 
供 了 基础 数据 。 
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Land use and land cover classifications of Horqin Sandy Land dune-meadow areas 


CAO Wenmei', LIU Tingxi, WANG Xixi, WANG Guanli', 
LI Dongfang', TONG Xing! 
(1. Water Conservancy and Civil Engineering College, Inner Mongolia Agricultural University, Inner Mongolia Water 


Resource Protection and Utilization Key Laboratory, Hohhot 010018, Inner Mongolia, China; 
2. Civil and Environmental Engineering, Old Dominion University, Norfolk 23529, Virginia, USA) 


Abstract: In this study, we determine land use/land cover (LULC) types using only optical remote sensing data 
in a dune-meadow area of Horqin Sandy Land in northeast China. We used 64 Sentinel-2 remote sensing images 
from 2018. An LULC decision tree recognition rule was established by combining image segmentation 
technology, vegetation phenology information, and habitat characteristics. The overall classification accuracy was 
0.91, and the Kappa coefficient was 0.89. The classification results show that most of the study region was 
dryland area, accounting for 33.79%, followed by shrub communities at 25.03%. High- diversity semi- shrub 
communities and arbor forests accounted for 14.54% and 10%, respectively, while low- diversity semi- shrub 
communities, meadowlands, and mobile sand lands account for about 5% each. The total proportion of other 
LULC types was less than 5%. The results show that this interpretation method better reflects the spatial 
distribution of the LULC while also providing growth and development data for different cover types. These data 
can be used to study material cycles and provide threshold references for historical LULC identification of 
Horqin Sandy Land. 

Keywords: remote sensing; land use/land cover; multi-resolution segmenation; decision-trees; classification; Sen- 
tinel 2-A/B; time series; Horqin Sandy Land 


